Magnetism and structure at a vacancy in graphene 
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The electronic structure, bonding and magnetism in graphene containing vacancies are studied 
using density-functional methods. The single-vacancy graphene ground state is spin polarized and 
structurally flat. The unpolarized state is non planar only for finite segments. Systems containing 
periodic arrays of vacancies displays magnetic transitions and metal-insulator transitions. 
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Introduction- Graphene is a two-dimensional (2D) 
sheet of carbon atoms forming a honey-comb lattice, 
with promising novel technologies and new physics P, 0- 
These include applications based on carbon nanotubes[3|, 
4] which are topologically equivalent to folded graphene 
ribbons. Graphene nanoribbons have delocalized states 
only in one dimension, with conductive and localized edge 
states. The localized states appear as "flat bands" near 
the Fermi energy, within tight-binding models which do 
not treat the lattice relaxation of the edge atoms 

When lattice relaxation is included via first-principles 
calculations, with the dangling bonds saturated with H 
atoms |8[, the electronic structure is found to depend 
specifically on the ribbon width and the relaxed edge 
structure. Band gaps are found to appear at the Fermi 
energy. The role of edges in graphene nanoribbons in 
producing magnetic ground states has also been studied 
using such methods [9(. 

At finite temperatures, entropy as well as the 2D 
structure of graphene favours disorder. These lead to 
a quasi-3D wavy conformation [Tol 11|. Many authors 



have also used vacancy or defect models [12j to include 
disorder effects in the graphene system. However, our 
calculations 13] for vacancies and also N-substituted de- 
fects in graphene showed that the energy costs of such 
defects are large enough to prevent natural occurences 
of such lattice defects, thus differing from metals or 
semiconductors. Hence in ref. [13) we concluded that 
the graphene samples used in initial quantum Hall-effect 
studies were good samples basically free of vacancies. 
However, vacancies and defects can be artificially intro- 
duced, and used to control the energy and spin states 
near the Fermi energy. The vacancies may be positioned 
lithographically to generate graphene containing a pe- 
riodic array of vacancies. Such arrays are metastable, 
and can be annealed to give various forms of tattered 
graphene. The tattering occurs because two second- 
neighbour vacancies are found to coalesce into a more 
stable divacancy, which form a sink for other second- 
neighbour vacancies. More distant vacancies simply dis- 
tort the cr-bonding skeleton and stabilize themselves. In- 
dividual vacancies carry a spin polarization, and provide 



a new type of qubit for quantum computational devices. 
The variations in spin polarization and band-gaps with 
the vacancy concentration set the stage for magnetic and 
metal-insulator transitions. Hence here we study an infi- 
nite graphene sheet with periodically distributed vacan- 
cies using a plane-wave density-function theory (DFT) 
approach based on the Vienna ab-initio simulation code 
(VASP 14]). We also study finite segments of graphene, 
with real-space DFT-calculations using Gaussian basis 
sets (Gaussian-98 code[l5j]). The finite graphene studied 
here have zig-zag edges terminated with hydrogen atoms. 

The finite graphcncs with a vacancy give a planar spin- 
polarized ground state, while the unpolarized state is 
slightly higher in energy (~ 200 meV) and show a non- 
planar, saddle-like structure (see Fig. QJ. 

The creation of the vacancy removes 4 electrons, and 
at the same time releases three sp 2 electrons (dangling 
bonds) and a 7r-electron. The sp 2 -type electrons are as- 
sociated with the atoms marked "a, b, c" in Fig. [1] The 
7r-electron is common to these three sites. The system 
distorts or spin polarizes to accomodate these dangling 
bonds and the free charges (4 per vacancy). The vacancy 
is surrounded by three pentagons. The tri-pentagonal 
structure has one pentagon smaller than the other two, 
enabling a singlet interaction between the atoms "b" and 
"c" which are closely positioned, compared to the pairs 
"a, b" and "a, c" of the other two pentagons. The atom 
"a" in the unpolarized structure (UPS), juts upwards 
(two views of the structure are given in Fig. [T]), while 
the atoms marked "b, c" dip downwards. This distortion 
disconnects the 7r-electron generated by the vacancy from 
the rest of the 2D electron network. The atom "a" of the 
UPS carries a lone-pair singlet of electrons, making the 
total structure unpolarized. The sheet around the va- 
cancy acquires a saddle-like distortion. The two bonds 
associated with 'X', and the atom 'a' defining the edge 
of the vacancy are shortened (~ 1.39A) compared to the 
normal 1.42 A bond length in graphene. 

The lower-energy (ground) state is the planar spin- 
polarized structure (PSPS). Two views of it are given 
in Fig. [TJ The stabilization energy for this system (44 
hexagons with a zig-zag edge saturated with H atoms), 




FIG. 2: The distribution of excess charge density (hatched 
areas) in the 2D plane of the graphene sheet around a vacancy, 
for the structure shown in Fig.[T] for the spin-unpolarized (a) 
and polarized (b) case. 




FIG. 1: (Color online)Relaxation of a finite sheet (44 
hexagons) of graphene containing a vacancy. The carbon 
atoms marked "a, b, c" were bound to the C atom removed to 
form the vacancy. The bent structure shows the non-planar 
graphene sheet where the four electrons released by the va- 
cancy are accomodated to form a spin-unpolarized structure. 
The ground state is planar and spin-polarized. Two views of 
it are shown. The atoms marked 'b', 'c' (in blue) are close 
enough in both structures to form a singlet interaction. The 
atom 'a' (in red) juts upwards in the unpolarized structure 
and carries an unpaired electron. The atoms 'X' and 'a' have 
two short bonds in the pentagons. 



compared to the unpolarized state is 217 meV, i.e., 
comparable to room-temperature energies. This result 
is obtained with a Gaussian-98 DFT calculation us- 
ing the Becke-Lee- Yang-Parr (B3LYP) exchange- corre- 
lation metafunctional[16| within a Gaussian (6-31G*, see 



Ref. [Uj for definitions) basis set. The short distance be- 
tween atoms "b, c" enables the two sp 2 dangling bonds to 
form a singlet. The atom "a" now carries one sp 2 -electron 
spin, unlike in the UPS where two (i.e., a lone-pair sin- 
glet) electrons are localized. Here again the bonds at "a" 
and "X" defining the edge of the vacancy are short. The 
7T electron released from the vacancy delocalizes into the 
2D network around the vacancy in the PSPS. In fig. [2] we 
show the distribution of excess charge projected on the 
x-y plane of the graphene sheet. 

On comparing the panels (a) and (b) of Fig. [5] it is 



clear that the charge depletion caused by the unpolar- 
ized structure is significantly larger than in the polarized 
structure of panel (b). In the PSPS, (b), the it electron 
released from the formation of the vacancy re-distributes 
itself among the three pentagons. This simplified dis- 
cussion of the 7r-electron charge at the vacancy has to 
be extended to allow the 7r-electron to spread in the full 
2D-network bounded by the edge-C atoms. The edge de- 
fines a a "zig-zag" termination (saturated with hydrogen 
atoms in this calculation). A plot of the Mullikan spin 
charge (n a — n^), obtained from a spin-polarized DFT 
calculation (using Gaussian-98) is shown in Fig.[3]for the 
structure of Fig. Q] The total spin of the structure is 
consistent with a triplet state made up of the sp 2 elec- 
tron on the carbon atom 'a' and the 7r-electron. It is 
also noted that some of the spin-polarization resides at 
the graphene-edge seen in the depth of the figure. This 
suggests that the finite-graphene segment studied here is 
not large enough for the vacancy to be isolated from the 
edge effects. The electronic density of states (DOS) for 
these are also shown in the bottom panel. The DOS of 
the unpolarized state was divided by 2 to compare with 
the polarized states, and a Lorentzian broadening of 0.3 
eV has been applied. The DOS shows substantial filling 
around Ep, due to the localized states at Ep created by 
the formation of the vacancy. 

The finite-graphenes necessarily limit the 7r- electrons 
to a localized region. The boundary conditions imposed 
on the wavefunctions do not allow true extended states. 
This difficulty is overcome in calculations where peri- 
odic boundary conditions in the x-y plane are used. We 
consider simulation cells with N a x Nt graphene unit 
cells, each cell containing two carbon atoms, with a lat- 
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FIG. 3: The top panel gives the spin density (n a — rip) where 
n a is the spin-dependent Mullikan charge on each C atom of 
the planar structure. The structure of the flat graphene piece 
is given in Fig.[Tj The bottom panel shows the Electronic den- 
sity of states (DOS) for the unpolarized system (bent struc- 
ture) and for the flat, polarized ground state. The Ef is set 
to zero and calculated with a Lorentzian broadening of 0.3 
eV. 



tice parameter an = 2.47A. and a C-C bond length 
of 1.42A in unperturbed graphene. Calculations are 
mostly for N a = Nf, for N a =2,3,4, and 6, correspond- 
ing to N = 2Nl of 8, 18, 32 and 72 carbon atom in 
the graphene simulation cell (SC). These calculations use 
the VASP[14| plane- wave code (clearly, no terminating- 
H atoms are needed) . Projected augmented- wave (PAW) 
pseudopotentials[14j have been used for Carbon. The C 
pseudopotential is well established and was used in sev- 
eral graphene- type calculations (e.g., Ref. [HI, 13, 18). 
Removal of one C atom from the test structures gives 
vacancy concentrations x v of 1/8, 1/18, 1/32, and 1/72. 
The finite-graphene piece studied using Gaussian-98 cal- 
culations contained 112 C atoms less one due to the va- 
cancy. Of these, 46 were edge atoms. Thus, the finite 
graphene is possibly analogous to an extended system 
with a vacancy concentration x v of ~ 1/66 if the edge ef- 
fect is negligible in the above finite-system. Similarly, the 
periodic simulation cell with 36 unit cells and 72 C atoms, 
x v = 1/72 is the most "isolated" vacancy studied here, 
the periodic image being six unit cells away. The systems 
with x v > 1/32 increasingly contain vacancy- vacancy in- 
teractions, since the periodic repetitions of the SC would 
place vacancies at 4, 3, and 2 units of an from each other. 

When vacancies are introduced into graphene, the 
large stresses are relieved by some neighbours (C-atoms 
near the vacancy) moving towards the vacancy, while 



others move away. The distortion persists to at least 
the third set of neighbours, and generates a bond-length 
distribution varying from about 1.37A to 1.45 A . If sev- 
eral vacancies are present in the system, the energy min- 
imization enables some vacancies to migrate and form 
divacancies. Thus the initial periodic structures can be 
"annealed" to give structures which are (at least) at local 
energy minima. Here we report results for systems with 
only one vacancy per simulation cell. The structure opti- 
mization involves (i) optimization of all atomic positions 
in a given simulation cell with N a x Nj, unit cells and 
optimization of the cell parameter an, (ii)repeating spin- 
density functional calculations till convergence of ener- 
gies, and negligible Hellman-Feynman forces (to less than 
0.001 cV/ Angstrom) are obtained. Details of integration 
grids, simulation cells etc., are similar to those discussed 
in previous work [ill ll3T |. Unlike in our results for finite 
graphene, no out-of-plane distortions were found in the 
periodic systems studied here, irrespective of the spin- 
polarization. The calculations for periodic systems also 
reproduce a tri-pentagonal vacancy, similar to the va- 
cancy in finite graphenes. Both systems are shown in 
Fig. [J 

In Fig. [4] we summarize our results for the magnetism 
and stabilization of graphene sheets containing periodic 
realizations of vacancies. The unpolarized planar system 
is the ground state for systems with vacancies in excess 
of x v ~ 0.06. At these higher vacancy levels, the atoms 
move sufficiently close to each other and singlet pairing 
becomes possible. At lower x v , the situation seems to be 
similar to the finite-graphene system where one of the sp 2 
electrons remains unpaired, localized at the atomic site 
'a'. The uncoupled 7r-electron becomes a part of the 2D 
electron fluid since the C-skeleton remains planar. The 
intriguing result obtained from the VASP calculation is 
that the net magnetization (i.e., n a — np ) accumulated 
within the simulation cell can takes a fractional value 
in the spin-polarized case. This result remains robust 
when tested by increasing the size of the fc-integration 
grids, and varying other technical parameters of the cal- 
culation. One possibility suggests some sort of spin and 
charge fractionation. The shortened pair of bonds at the 
atom 'a', and at the atom 'X' suggests that additional 
electron localization occurs around them. The 7r-electron 
released in forming the vacancy remains near the vacancy, 
shared among the three pentagons defined by the atoms 
'a', 'b' and 'c' of the vacancy (Fig. [1]). However, a 1/3 
spin-charge seems to be available for alignment with (or 
against) the spin already localized on the atom 'a'. The 
calculations attempt to examine such possibilities using 
explicit spin-density functional calculations which reli- 
ably predicts spin ground states. It is found that the 
spin localizes around the vacancy and aligns ferromagnet- 
ically for x v < 1/20, (given the maximum spin of 1.333) 
and antiferromagnetically for 1/20 < x v < 1/16, i.e., a 
magnetization of 0.66). Thus the low vacancy regime 
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FIG. 4: (Color online)The stabilization energy E po i — E unp , 
and the spin polarization of the graphene vacancy, with one 
vacancy per simulation cell (containing N lattice sites), and 
periodically repeated. The polarized and unpolarized struc- 
tures are planar. The unpolarized planar ground state is sta- 
ble for x v >~ 0.06. 

gives a net spin polarization n z = (n a — rip) such that 
1.333 < n z < 0.666, while the intermediate regime shows 
0.666 < n z < 0. The high vacancy regime is unpolarized. 
The present calculation cannot sharply pin down the ac- 
tual transitions or rigorously establish these suggestions. 

The electronic DOS indicated in Fig. [3] for finite 
graphene shows significant density of states near Ep- The 
situation is similar for infinite graphene with periodic va- 
cancies. The high vacancy (unpolarized) systems have 
an energy gap. However, depending on the location of 
Ep within the band structure, even within a rigid band 
model we see that the system may or may not be a metal, 
depending on the vacancy concentration x v . This issue 
needs a fuller discussion in a separate study. 

The fore-going considerations suggest that the inter- 
action of the 2D electrons containing the 7r-electron re- 
leased from the vacancy, with the localized sp 2 -electron 
spin on the atom 'a' has features of the Kondo problem, 
as well as fractional-charge/spin physics. Investigation 
of such matters would be best carried out within a more 
analytical approach, complemented by the results of com- 
putational approaches similar to the present study. 

In conclusion, we have demonstrated, using finite as 
well as periodic structure calculations, that the ground 
state of a graphene sheet (held in free space at zero tem- 
perature) containing an isolated vacancy is spin polar- 
ized. When the concentration of vacancies (realized as a 
periodic array) exceeds ~ 0.06, the unpolarized structure 



becomes the ground state. All calculations using periodic 
boundary conditions lead to planar structures. Only the 
finite-graphene pieces which are unpolarized become dis- 
torted. 
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